Quantum wire junctions breaking time reversal invariance 
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We explore the possibility to break time reversal invariance at the junction of quantum wires. The 
universal features in the bulk of the wires are described by the anyon Luttinger liquid. A simple 
necessary and sufficient condition for the breaking of time reversal invariance is formulated in terms 
of the scattering matrix at the junction. The phase diagram of a junction with generic number of 
wires is investigated in this framework. We give an explicit classification of those critical points 
which can be reached by bosonization and study the interplay between their stability and symmetry 
content. 



I. INTRODUCTION 



Time reversal symmetry is a fascinating subject. 
In this paper we investigate the behavior of junc- 
tions of quantum wires under time reversal transforma- 
tions. Quantum wire networks with junctions, which 
attract recently much attention 1-29 , are essentially one- 
dimensional systems whose transport properties are af- 
fected by quantum effects. The universal features in the 
bulk are captured by the Luttinger liquid theory 30 . The 
junctions represent in this context a kind of quantum im- 
purities (defects) , where both reflection and transmission 
can take place. This fact gives origin of a complicated 
phase diagram, which has not been yet fully understood 
for general boundary conditions at the junctions, formu- 
lated in terms of the basic fermion fields. Focussing on 
the case of one junction, we discovered 19 in the frame- 
work of bosonization a large class of boundary conditions, 
which preserve the exact solvability of the Tomonaga- 
Luttinger (TL) model describing the Luttinger liquid in 
the bulk. At criticality these boundary conditions simply 
express the splitting of the electric current in the junction 
and are therefore quadratic in the fermion fields. We clas- 
sified and studied in this setting all critical points which 
respect time reversal invariance. In this paper we ex- 
tend our framework in order to cover also that part of 
the phase diagram, where the time reversal symmetry is 
broken. Recalling that the Tomonaga-Luttinger dynam- 
ics preserves time reversal invariance, the breaking can 
take place only at the junctions. In principle such kind of 
junctions can be realized 9,10 ' 15,24 by means of an external 
magnetic field and are therefore of practical interest. 

The previous theoretical investigations of the stability 
of the critical points and their behavior under time rever- 
sal have been mostly focussed on junctions with n = 3, 4 
wires. Applying the framework developed in Refs. 19-21, 
we face below these problems for generic n. 

The paper is organized as follows. In the next section 
we define the bulk dynamics and boundary conditions at 
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FIG. 1: A star graph T with n edges modelling the junction 
of quantum wires. 



the junction. Using bosonization, we recall 21 in section 
III the exact (anyon) solution of the model. In section IV 
we derive the current-current correlation function and ex- 
tract the necessary and sufficient condition for the break- 
ing of time reversal. We discuss here also the Kirchhoff 's 
rules relative to the U(l) ® U(l) symmetry of the model. 
In section V we consider the conductance and describe 
the impact of time reversal breaking on it. The classifi- 
cation and parametrization of the critical points is done 
in section VI. In section VII we study the phase diagram, 
concentrating mainly on the symmetry content and sta- 
bility of the fixed points. Section VIII is devoted to our 
conclusions. Some technical details are collected in the 
appendices. 



II. BULK DYNAMICS, SYMMETRIES AND 
BOUNDARY CONDITIONS 



The quantum wire junction is modeled by a star graph 
r of the form shown in FIG. 1. The edges Ei are half- lines 
and each point P in the bulk of T is uniquely determined 
by its coordinates (x,i), where x > is the distance 
to the vertex V and i = 1, ...,n labels the edge. T \ V 
represents the bulk of the graph. The bulk dynamics is 
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governed by the TL Lagrangian density 

C = hj)l{dt - v F d x )ip! + iip 2 (d t + v F d x )ip 2 

- g+tyWi + ^2? - g-ffl1>i - ^) 2 ■ (1) 

Here {ip a (t, x,i) : a = 1, 2} are complex fields, vf is the 
Fermi velocity and j±el are the coupling constants 31 . 

The bulk theory has an obvious £/(l)®i/(l) symmetry. 
In fact, the Lagrangian density (1) is left invariant by the 
two independent phase transformations (s,ssK) 

1>a - e is Va , C - e" is V* , (2) 

^ - e-^" 1 )" 8 ^ , C - e^ 1 )"^ , (3) 
implying the current conservation laws 

d t p± (t, x,i)- v F d x j± (t, x,i)=0, (4) 

where 

p±(t,x,i) = (V'iV'i ± $2^2) {t,x,i) (5) 

are the charge densities and 

j±(t,x,i) = p T (t,x,i) (6) 

are relative currents. We adopt below also the chiral 
combinations 

j R (t,x,i) = ^(C-i- +C+j+)(t,x,i) , (7) 

j L (t,x,i) = ^(C-i- - C+i+)(*.a ; ,*) . (8) 

where the real parameters (±, determined later on, arc 
such that jl and j R represent the particle excitations 
moving towards and away of the vertex respectively. In- 
terpreting the vertex as a defect, characterized by some 
scattering matrix, the currents jl and ]r describe there- 
fore the incoming and outgoing flows. 

The bulk theory is invariant also under the time rever- 
sal operation 

T^ x (t,x,i)T* = 4> 2 {-t,x,i), (9) 
Ti> 2 (t,x,i)T* = ^{-t,x,i), (10) 

where T is an anti-unitary operator. As already men- 
tioned, our main goal below will be to investigate the 
impact of the vertex V and the related boundary condi- 
tions on time reversal. 

The TL model (1) is exactly solvable on the line K, 
but much care is needed on the graph T, where some 
boundary conditions must be imposed at the vertex V. 
Keeping in mind that the quartic bulk interactions in (1) 
can be solved exactly via bosonization 30 , it will be ob- 
viously convenient to formulate the boundary conditions 
directly in bosonic terms. In this spirit and according 
our previous comments on the chiral currents, it is quite 
natural to require that at a critical point 

n 

j L (t,0,i) = J2 s ikJR(t,0,k), VteR. (11) 

k=l 



For n — 2 this boundary condition has been first pro- 
posed and explored in Ref. 2. Because of scale invariancc 
at criticality, § is a constant (momentum independent) 
unitary scattering matrix, 

SS*=I. (12) 

Since the chiral currents (7,8) are Hermitian fields, one 
requires also that S has real entries, 

S = S. (13) 

Eqs. (12, 13) imply that S is any element of the orthog- 
onal group 0(n). It has been shown in Refs. 19-21 that, 
in spite of the fact that the boundary condition (11) is 
quadratic in the fields tp a , it preserves the exact solvabil- 
ity of the TL model on the graph T. It is worth men- 
tioning that this is not the case with the linear boundary 
conditions in ip a , which might look at first sight simpler. 

Applying the time reversal operation (9,10) to (11) one 
infers 

n 

j L (t,0,i) = J2§ik3R(t,0,k), VteM, (14) 
fc=i 

where the apex t stands for transposition. Comparing 
(11) and (14) we conclude that symmetric scattering ma- 
trices 

S = S* (15) 

define boundary conditions which respect the time rever- 
sal invariance. This is the case we investigated previously 
in Refs. 19-21. On the other hand, for 

§ ^ S* (16) 

one expects breaking of time reversal. We demonstrate in 
Section IV that this is indeed the case, using the explicit 
form of the current-current correlation function. 

We conclude at this point the concise description of 
the bulk dynamics, symmetries and boundary conditions 
of our model and briefly describe in the next section the 
solution. 



III. SOLUTION OF THE TL MODEL ON A 
STAR GRAPH 

We look below for the solution ip a of the TL model 
which satisfies the boundary condition (11) and obeys 
the anyon exchange relations 

i()*(t,xi,i)i) a (t,X2,i) = 

e (-i)-i^(x 12 )^ a(t) X2 ,w a (t, Xl ,i). (17) 

Here e(x) is the sign function, x 12 = X\ — x 2 and n e M is 
the so called statistical parameter, which equals an even 
and an odd integer for bosons and fermions respectively. 
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Other values of k give rise to Abelian anyon statistics 
"interpolating" between bosons and fermions. 

The solution on T can be expressed in terms of the 
chiral scalar fields {fi^z{C> ■ Z = L,R; i = l,...,n}, 
which are not independent as on the line K, but respect 
the constraints 



mAZ) > 



(18) 



keeping track of the boundary conditions (11). The ex- 
plicit construction and a summary of the main features 
of tpi t z are given in appendix A. A key point is the non- 
trivial one-body scattering matrix 32 



S(k) = 0(-k)S + 0(fc)S* , 



(19) 



where 9 is the Heaviside step function. We stress that the 
peculiar /c-dependence of S(k) respects scale invariance. 

Let us summarize now the basic features of the solu- 
tion of the TL model with boundary conditions (11). We 
do this essentially for two reasons. First of all the field 
<f associated with the S'-matrix (19) behaves quite differ- 
ently (see appendix A) from its counterpart in Ref. 21. 
Second, because we would like to keep the present paper 
self-contained. 

Following the standard bosonization procedure 30 , we 
set 



ipi(t, x, i) — Zi : e 1 
ip 2 (t,x,i) = z l : e 



spii [a If i , fi .( V t - X ) + T tp i , L ( V t + X ) ] 

\/Tr[Tcpi tR (vt-x)+cripi iL (vt+x)] 



:,(20) 
= ,(21) 



where : • • • : denotes the normal product relative to the 
creation and annihilation operators of the fields <fi,z, 
namely the generators of the algebra (A2) . The explicit 
form of the normalization constants z% (including the 
Klein factors) is reported in appendix A as well. Finally, 
a, t and v are three real parameters to be determined 
in terms of coupling constants g± and the statistical pa- 
rameter k. For this purpose we can assume without loss 
of generality that 

cr>0, a^±T (22) 
and introduce for convenience the variables 

C±=r±<7. (23) 
Plugging (20,21) in (17) one gets 

C+C-=k- (24) 

Moreover, using standard short-distance expansion for 
the charge densities, one obtains 

p±(t,x,i) = 

— ^— [(d<p itR )(vt -x)± {d<pi, L )(vt + x)} , (25) 



the normalization being fixed 21 by the U(l) (8) U(l)-Waxd 
identities. Inserting (20,21,25) in the quantum equations 
of motion 

i[dt + {-l) a v F d x ]ip a (t,x,i) = 

2[.9+ : P+(t,x,i)ip a : {t,x,i) 

-(-l) a 9- :p-(t,x,i)i> a : {t,x,i)}, (26) 



one finds 



v f k + -g + 

7T 

2 

vpK H — g_ 

7T 



(27) 
(28) 



Eqs. (24,27,28) provide a system for determining v and 
C± (equivalently a and r) in terms of vp and g±. The 
solution is 



\ttkvf + 2g- 



(29) 



v = V(kkv f + 2gJ)(-KKV F + 2g + ) 



7T K 



where the positive roots are taken in the right hand 
side. The relations (29) and (30) represent the anyonic 
generalization 21 of the well known result for canonical 
fermions (n = 1) in the TL model 33 . The conditions 
2g± > — ttkv f ensure that a, r and v are real and finite. 

Finally, in the bosonic variables U(l) <8> [/(l)-currents 
j± take the form 



3±(t,x,i) 



2^kv f Q 



[(d<pi, R )(vt -x)t {d<Pi, L ){vt + x)] (31) 



and satisfy (4) by construction. Using (18) and (31) one 
immediately verifies that the above solution of the TL 
model on V indeed satisfies the boundary condition (11). 



IV. SYMMETRY CONTENT 
A. Time reversal 

The simplest way to investigate the behavior of the 
above solution under time reversal is to derive the two- 
point correlation functions of the currents j± , defined by 
(31). Using (A10-A12) one obtains 



(j+(t 1 ,X 1 ,i 1 )j+(t 2 ,X2,i2)) 



(2ir(+v F ) 2 



[S ili2 T> 2 (vt 12 - x 12 ) + 5 lll2 V 2 (vt 12 + x 12 ) 



- S lll2 V 2 (vt 12 + x 12 ) - §l it V 2 (vt 12 - x 12 )\ , (32) 



where 



(33) 
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and ti2 = ti - t 2 , x 12 — x\ - x 2 and x X2 = x\ + x 2 . 

Let us assume for a moment that time reversal is an 
exact symmetry, or equivalently that T leaves invariant 
the vacuum state fi. Then, using the anti-unitarity of T, 
one finds that 

(j + (t 1 ,x 1 ,i 1 )j + (t 2 ,x 2 ,i 2 )) = 



, x x , h)j + (-t 2 , x 2 ,i 2 )} (34) 

holds. Combining (32) with (34) one deduces that 

TQ = Q § = S* , (35) 

showing that the TL model on T is invariant under time 
reversal if and only if § is symmetric. Otherwise, time 
reversal is broken, i.e. 



(36) 



which confirms the conjecture made after equation (16) 
in the introduction. In particular, time reversal is always 
exact for n = 1 . For this reason we focus in what follows 
on the case n > 2. 



B. U(l) (g> (7(l)-symmetry 

Continuous symmetries on graphs are governed 19 ' 20 by 
the associated Kirchhoff's rules. Concerning the U(l) <g> 
[/(l)-symmetry, using the current conservation (4), one 
gets for the corresponding charges 

n „oo n 

d t Q± = dt'Y] I dxp ± (t, x, i) = v F y~] j± {t, 0, i) . 
i=i Jo i=1 

(37) 

Inserting here (31) and taking into account the boundary 
conditions (18), one finds 



n 



(38) 



V. CONDUCTANCE 

A simple physical observable, which is sensitive to the 
breaking of time reversal, is the conductance tensor 
of the Luttingcr liquid on T. In order to compute this 
tensor, one couples the theory to an external potential 
A x (t,i) by means of the substitution 



d x i — ► d x + iA x (t, 



(40) 



in cq. (1). The resulting Hamiltonian is time dependent 
and the conductance is the coefficient in the linear term 
of the expansion of the expectation value (J x (t,0,i))A x 
in terms of A x . For deriving G one can apply therefore 
linear response theory, which leads to 19,20 



2ttv f C 



-(I- 



(41) 



Using the condition (39), which ensures the conservation 
of the electric charge Q + and S € 0(n), one gets the 
Kirchhoff's rule for the conductance tensor 

n n 

^Gy= 5^0^ = 0, Vj = l,...,n. (42) 

i=l i=l 

If a voltage Vi is applied to the edge E i} the current Ij 
flowing in Ej is 



Ij — ^2 <&jiVi ■ 



(43) 



Combining (16) with (41) we conclude that the breaking 
of time reversal (16) implies the asymmetry 



(44) 



a feature which has been previously observed in Rcfs. 
7,10,22. The property (44) provides an attractive ex- 
perimental signature. Indeed, consider for instance the 
following two configurations with i ^ j. Apply first the 
voltage V to the edge Ei, setting to the voltages in all 



other edges, and measure the current Ij 



iV. Rc- 



From this result wc infer that 19 , 



peat the same operation, applying now V to the edge Ej 
and measuring the current Jj = G%jV. If h/Ij ^ 1, the 
system breaks time reversal. 



Q± — conserved ^ Sjj = ±1 Vj = 1, n . (39) 



VI. CRITICAL POINTS 



Recalling that Q+ is the electric charge and Q_ the 
hclicity of the Luttinger excitations, we see that only 
one of these quantum numbers is preserved for a generic 
junction 34 . 

It is worth mentioning that for junctions with n = 2 
wires the conservation of Q + or Q_ protects the time 
reversal symmetry. Thus, junctions of three wires (T- 
junctions and Y-junctions) represent the minimal setting 
for breaking time reversal in systems preserving the elec- 
tric charge Q + or the helicity Q_ of the Luttinger liq- 
uid. Notice also that the conservation of Q + excludes 
the Dirichlct fixed point S = — I. 



A. Classification 

As already mentioned, § € 0(n). There exists there- 
fore an orthogonal matrix O, such that 

o 



±i 



(45) 



■5 



Here r 4 are q rotation matrices 

f cos 6i — sin#A . , . ,.„. 

ri ={ S inel cob* J" (46) 

Let us denote by p± the number of eigenvalues ±1 of 
(46). Then q = ^(n — p + — p~) and the critical points 
are classified by the set (p + ,p-, 8\, ...,8 q ). From (36) we 
conclude that the time reversal symmetry is broken if 
and only if 8k ^ — 7r, for some k = l,...,q. The angles 
6k thus codify the breaking of time reversal. 

B. Parametrization 

§ can be any element of 0(n), but in the physical appli- 
cations one is mostly interested in boundary conditions 
which preserve the electric charge Q+. In this case one 
infers from (39) that 

§v = v, v=-^=(l,l,...,l), (47) 

I 



i.e. § leaves invariant the vector v. Let R be the orthog- 
onal matrix (see equation (B2) in appendix B), which 
rotates the vector (0, 0, 0, 1) in v. Than § admits the 
representation 



S = R ^ f) R* (48) 



with S' € 0(n — 1). Therefore the boundary conditions 
which respect the [/(l)-symmetry of the TL model are 
parametrized by the group 0(n— 1). The two connected 
components of 0(n — 1) give origin of two continuous 
families of critical points. Each family depends on (n — 
l)(n — 2)/2 parameters, which are the angular variables 
parametrizing the elements of 0(n — 1). 

Let us illustrate this simple general structure for n = 3. 
The two families of critical points depend in this case on 
one angle d G [— 7r,7r) and read 



1 + 2 cos j? 1 — cos tf + V3 sin i? 1 — cos — y/3 sin z9 

1 — cos "d — y/3 sin i9 l + 2cos , i9 1 — cos <? + V3 sinz9 

1 - cos + Vising 1 - costf - V3sintf l + 2cos$ 



(49) 



, , l-2cost? 1 + cos?9 - Vising 1 + cos ft + Vising 

§ (2) (^)=o l + cos^-V3sini9 1 + cos i? + V3 sin l-2costf 

' 1 + cos -d + Vising l-2cosi? 1 + costf - \/3sinz9 



(50) 



confirming the recent results of Ref. 26. The point 



-1 2 



(-7T) 



-1 2 



(51) 



In appendix B we report an explicit parametrization 
of the n x n critical S-matrices. The case when the (7(1)- 
symmetry is preserved can be treated analogously 20 . 



has been discovered by Griffith 35 more than five decades 
ago in his pioneering work on graph models in quantum 
chemistry. According to (41), in this case the conduc- 
tance of the Luttinger liquid is enhanced with respect 
to the line, which has been associated 3 with the phe- 
nomenon of Andreev reflection. The Neumann point 



§ (1) (0) 




(52) 



describes instead an ideal isolator because G = 0. To 
our knowledge the whole family §( 2 )($) was derived 36 
first in Ref. 19 and, together with the points (51, 52), 
preserves time reversal invariance. Finally, the matrices 
S^{^) with d 7^ — 7r,0 give all critical points which vi- 
olate time reversal symmetry for n — 3. In a different 
parametrization 37 they appeared in Refs. 10,22. 



VII. PHASE DIAGRAM: BOUNDARY 
DIMENSIONS AND STABILITY OF CRITICAL 
POINTS 

The boundary dimensions of the solution ifj a capture 
the impact of the junction at criticality and can be ex- 
tracted from the two-point functions 

(ipl(t 1 ,x 1 ,i 1 )ip 1 (t 2 ,x 2 ,i 2 )) = z h z i2 

[V(vt 12 - x X2 )\ aH ^ [V(vt 12 + x 12 )Y 25 ^ 
[V(vt 12 - x 12 )r s ^2 [V(vt 12 + x 12 )r s ^ , (53) 



and 



(ip 2 (ti,xi,i 1 )ip 2 (t 2 ,x 2 ,i 2 )) = (53) with a <-> r . (54) 
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Performing the scaling transformation 

t *-^> gt , x^> qx , g > . 
in (53, 54), one obtains 



(55) 



(tpaigh, gxi,ii)tp a (gt 2 , gx 2 ,i2)) = 

g - D n*2 (ip*(ti,x 1 ,i 1 )ip a {t 2 ,x 2 ,i 2 )} , (56) 



where 



D = (a 2 +T 2 )l n +aT(S + § t ) 



(57) 



The eigenvalues di of the matrix D/2 are the scaling di- 
mensions. If time reversal is broken (S 7^ §'), some of the 
eigenvalues of S are necessarily complex. Notice however 
that the eigenvalues of the combination S + §* are all real 
and 



di = i((T 2 + T 2 ) + (TTS, , 



(58) 



s being the n-vector 



s = (cos 61, cos 6\, cos 6 2 , cos 2 , cos6 2q , ±1, ...,±1) . 

(59) 

In appendix C we prove that the mixing between ipi and 
i/j 2 produces vanishing additional eigenvalues and there- 
fore does not affect the spectrum (58,59). 

Recalling that the scaling dimension on the line is 38 



d^ = l -{a 2 +r 2 ), 
one deduces from (58) the boundary dimensions 



(boundary) = 



e + -e 



Si . 



(60) 



(61) 



which control 38,39 the stability of the critical points. The 
direction i at a critical point § of the phase diagram is 
stable (unstable) if di > (di < 0). We call the point § 
completely stable if all relative directions are stable. Us- 
ing (27,28), the boundary dimension di can be rewritten 
in our case in the form 



,(boundary) 



2^+ 



-)8i 



v>0. 



(62) 



where v > is given by (30). It is natural to consider 
at this point the two regimes of repulsive (g + > g_) and 
attractive (g + < g_) anyonic interactions. From (62) 
one concludes that in the repulsive case the direction i is 
stable if Si > 0. Vice versa, in the attractive case stability 
requires < 0. 

The direction i in the phase diagram is called flat if 
di = 0. This happens for cos#j = and/or g + = g_. The 
last case is very special: there is no interaction between 
and tp2 (see eq. (1)), all boundary dimensions vanish 
and all directions are flat. 

It is worth stressing that the above considerations con- 
cern the phase diagram of the system without symmetry 



constraints. According to Section IV however, the Kirch- 
hoff's rules controlling the symmetry content of the TL 
model on T, impose such constraints. If one requires for 
instance [/(l)-symmetry, the condition (39) implies that 
Si = 1 in at least one direction. Therefore, for attractive 
interactions with U(l) symmetry there are no completely 
stable points. The same conclusion holds in the repulsive 
case with U(l) symmetry. 

As already mentioned, imposing time reversal symme- 
try implies that 9i — —it or 6i = for all i = 1, ...,q, 
which severely restricts the phase diagram. In partic- 
ular, the only completely stable fixed points arc § = I 
(for g + > and S = — I (for g + < <?_), corresponding 
to Neumann and Dirichlet boundary conditions respec- 
tively. Allowing for breaking of time reversal leads to a 
richer phase diagram, which admits whole families of non 
trivial (S 7^ ±1) completely stable critical points. 

Let us consider for illustration the phase diagram for 
n = 3 (Y-junction). We have shown above that all critical 
points, respecting the electric charge conservation, are 
given by (49, 50). The corresponding eigenvalues are 



(cos 0, cos 0,1), s (2) = (1,1,-1) 



(63) 



showing that the family §^(9), which preserves time 
reversal symmetry, does not contain completely stable 
points. The time reversal breaking family SW(#) con- 
tains instead the non trivial completely stable fixed 
points with cos 8 > in the repulsive regime and cos 9 < 
in the attractive one. In this sense complete stability 
is favored by time reversal breaking. 

We stress in conclusion that the above algorithm can 
be applied for analyzing the stability of the critical 
points under perturbation with any composite operator 
involving the basic fields ip a (t,x,i). Some examples of 
quadratic operators are considered in appendix D. 



VIII. CONCLUSIONS 

We investigated above the behavior under time rever- 
sal of a Luttinger junction with any number of edges 
n and satisfying the boundary conditions (11). As ex- 
pected, time reversal invariance can be broken by bound- 
ary effects, in spite of the that fact that the bulk theory 
preserves this symmetry. The following two exceptions 
are worth mentioning. Time reversal symmetry is al- 
ways preserved for n = 1 . The same conclusion holds for 
n = 2, provided that the electric charge Q + is conserved. 

The results of this paper give a global view on the 
phase diagram of the system with boundary conditions 
(11) and the framework allows to investigate both the 
symmetry content and the stability of the critical points. 
It turns out that the phase diagram has two connected 
components, corresponding to those of the group 0(n) 
and therefore depending on n(n— 1)/2 parameters, which 
describe irrelevant boundary couplings. In this classifica- 
tion the critical points, which respect the electric charge 
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FIG. 2: A graph with 3 external and 3 internal edges. 

conservation, form a 0(n — l)-subfamily. A simple cri- 
terion (16) allows to distinguish the points which vio- 
late time reversal invariance from those which preserve 
it. The stability of the critical points is controlled by 
the relative boundary dimensions. For generic n we de- 
rived these dimensions in explicit form (62), establishing 
their dependence on the boundary conditions and the 
bulk couplings. The analysis of the critical points, which 
are stable in all directions of the phase diagram, reveals 
that except of the Neumann point § = I for repulsive 
interactions and the Dirichlet point § = — I in the at- 
tractive case, all other completely stable points violate 
time reversal invariance. 

As already mentioned in the introduction, the simplest 
realization of devises, violating time reversal invariance, 
uses 9,10 ' 15,24 magnetic fields. An example, which fre- 
quently appears in the literature 10,22 , is the configuration 
shown in Fig. 2. One has three external half lines and a 
ring composed of three compact internal edges and three 
junctions. A magnetic flux (f> is crossing the ring. The 
complete field theory analysis of the Luttinger liquid on 
a graph with this geometry is very complicated problem, 
which is beyond the scope of the present paper. One ap- 
proximate way to face the problem could be to use the 
star product approach 40 or the "gluing" technique 41 ' 42 
for deriving the 3x3 scattering matrix relative to the 
external edges. Although a bit complicated 40,42 ' 50 , this 
S-matrix can be used afterwards for developing a simpli- 
fied model with one effective junction. Clearly, such an 
approach does not provide the conductance of the inter- 
nal edges Ii. 

The generalization of the results of this paper to off- 
critical junctions represents also a challenging open prob- 
lem. The study of the rich spectrum 29 of effects away 
of equilibrium is essential in this respect. Another in- 
teresting subject is the study of networks with several 
junctions. We are currently investigating these issues. 

Acknowledgments 

We thank Pasquale Calabrese for his interest in this 
work and for fruitful discussions. Correspondence with 
Ines Safi is also kindly acknowledged. The research of 
B.B. has been supported in part by the NSF grant PHY- 



0757868. 



APPENDIX A: CHIRAL FIELDS ON T 

The chiral scalar fields 

<p ilR (0 = / —7= K*(fc)e ifc « + «,(fc)e- ifc «] , 

JO 7TV2fc 

^ L (0 = / —7= [a*(-k)c^ + ai (-k)c-^] , 
Jo 7tv2k 

(Al) 

are the building blocks of the solution (20,21)). On T the 
generators {a,(fc), a*(fc)} obey the following deformation 

[ai(k) , aj(p)] = [a*(k) , a*(p)} = , 
[oi(fc) , a»] = 2ir[S(k - p)5 ij + 5(k + p)S ij (k)\ , (A2) 

of the standard canonical commutation relations. Here 
S(k) is the one-body scattering matrix defined by (19). 
Besides (A2), we impose also the constraints 

n 

a t {k) = J2 S ij( k ) a j(- k )> ( A3 ) 

n 

a*(k) = £ a *(-fc)%(-fc), (A4) 
i=i 

which are consistent, because S(k)S(—k) = I, and im- 
ply (18). Equations (A2-A4) define a special reflection- 
transmission algebra A, which has been introduced in a 
more general form in the study 43-46 of point-like defects 
in integrable systems. Notice that although fc-dependent, 
S(k) is scale invariant. 

Time reversal is realized in the algebra A by means of 

Ta t (k)T* = -ai(-k) , Ta*(k)T* = -a*(-k) . (A5) 

In fact, (A5) imply 

Tpi, R (t - x)T* = -<p itL (- t + x), (A6) 
T<pi, L (t + x)T* = -<p itR (-t- x ), (A7) 

which implement in turn the time reversal transformation 
(9,10) on the solution (20,21). 

For the construction of correlation functions we adopt 
the Fock representation of A. We denote by and (• , •) 
the Fock vacuum state and the scalar product, using for 
the vacuum expectation values of the operators Ok the 
short notation 

(n,o 1 ---o n n) = (o 1 ---o m ). (as) 

Since ai(k)Q = 0, the basic correlators are 

(ai(p)a*(q)) = 2n [% S(p - q) + S y (p) 5(p + q)\ , 
(a* i (p)a j (q))=0, 

(A9) 
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which imply 



(Vii,fl(^i)Vi2,ii(6)) = 

(<Pil,i(£l)<Pi2 = S»i» 2 «(/^12) , 

(^!,i?(6)^2,i(6)) = S- li2 u(^i2), 



where £ i2 = 6 - £2, 



u(0 = hi(i£ + e) , e > , 

7T 



(A10) 

(All) 
(A12) 



(A13) 



and n > is an infrared mass parameter 47 . The normal- 
ization constants Zi which occur in (20, 21) depend on /i 
in the following way 



(2^)- 1 /V [(ct2+t2)+2cttS " 1/2 ^ , 



(A14) 



where % are the anyon Klein factors needed to ensure 
the correct anyon exchange relations on different edges 



J 



of the graph T. A simple representation is 
rn = : e^+O : , 



(A15) 



where {a*, cti : i = 1, ...,n} generate the auxiliary alge- 
bra 



[ttj , aj] = [a* , a*] = , [a* , a*] = i-ey , (A16) 



with e. 



T for i < j , cu = and = 1 for i > j. 



It is worth stressing that there is an action principle 
behind the whole structure (A1-A13). The action can be 
written in terms of the combinations 

<Pi{t,x) = ^[(p itR (t - x) + ifi, L (t + x)] , (A17) 
&{t,x) = ^[(Pi,R(t-x)-ipi, L (t + x)] (A18) 

and the auxiliary fields {Xi(t, x), Aj(t, a;)} as follows. The 
bulk and boundary actions action are 



'bulk 



d* 

00 Jo 



+ 00 



dx 



i=l 



[A;(<9 X <^ + d t ipi) + Xi(d t ipi + d x tpi) (t, x) , 



(A19) 



^boundary 



^ n /■+00 



i=l 



n r+00 

- Y, d< t^( § + S *)«& - y*( s + s *)«Vj - 2 ^( § - s %£?] (t, 0) : 



(A20) 



respectively. The total action S = <Sbuik+<Sboundary 1S non 
degenerate, A and A being the conjugate momenta of ip 
and !p respectively. In agreement with (36), the only term 
breaking time reversal invariance is the term proportional 
to S — S in (A20). The bulk variation involves only <Sbuik- 
Varying with respect to A and A, one gets the duality 
relations 

d t <p(t,x,i)=-d x <p(t,x,i), (A21) 
d x <p(t,x,i)=-d t <p(t,x,i). (A22) 

The bulk variation with respect to ip and <p gives analo- 
gous relations between A and A. The boundary variation 
involves both Sbuik and ^boundary and, as easily verified, 
generates the boundary condition (18). 

A final comment concerns an interesting interplay be- 
tween locality and time reversal symmetry on T. A 
standard computation shows that at space-like separated 



points t\ 2 < x\ 2 



[(p(t 1 ,x l7 i) , tp(t 2 ,x 2 ,j)] = 

- [<p(h,xi,i),<p(t2,X2,j)] = i (S*-S) tf , (A23) 



implying that the time reversal breaking on T is accompa- 
nied by violation of locality of ip and (p. 48 One can easily 
check however that this violation does not affect the lo- 
cality of the currents j± , which belong to the observables 
of the theory. 

More about quantum field theory on graphs (also away 
of criticality) can be found in Rcfs. 18-20, 28 and 41- 
50, where some basic elements 51 ' 52 of the spectral theory 
of differential operators on graphs ("quantum graphs"), 
have been used. 
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APPENDIX B: CRITICAL S-MATRICES FOR 
GENERIC n 

First of all, the matrix R which rotates the vector 
(0, 0, 0, 1) in v can be taken in the form 

'0 ifi< j = l,2,..,n-l, 

if i > j = l,...,n- 1, 

if i=j = 1, ...,n- 1 , 
if i = 1, ...,n, , j = n, 

(Bl) 

As well known, the matrix §' G 0(n — 1) can be 
parametrized in terms of the (n — l)(n — 2)/2 rotation 
matrices {^(^i?) : i,j = 1, ...,n— 1, i < j} each of them 
rotating at the angle in the ij-plane. If det(S') = 1 
one has 



Rij — < 



y/ (n-j)(n -j + l) 

j n—i 
V n-i+1 



S' 



n 

Vi=n-2 



' i,n— 1 




n-3 



(B2) 

The only delicate point is the domain of the generalized 
Euler angles $ij , which turns out to be 53 

n ._/[-*■,*■) forj = i + l, 

\[-7r/2,7r/2] forj>i + l. { 6) 

Finally, in the case dct(S') = — 1 one can simply multiply 
the right hand side of (B3) by the matrix r, which reflects 
for instance along the first axis. 



APPENDIX C: THE Vi - 1P2 MIXING 

The mixing between ipi and V2 is described by the 
two-point functions 

{i>i(h,xi, 11)^2^2, X2, 12)) = Zi ± z i2 

[D(vh2 - x 12 )f rS ^ [V(vt 12 + x 12 )Y tS ^ 

[D(vt 12 - iu)]^ [D(vt 12 + xi2)V 2S ^ , (CI) 

and 

{i>2(tl,X 1 ,i 1 )i) 1 (t2,X2,i2)} = (CI) with (7 <-> t . 

(C2) 

Combining (C1,C2) with cqs. (53,54), one finds that 
under a scaling transformation (55) a generic two-point 
function transforms according to 

(d (eh, QXi,h)ip a2 (gt2, 9x2,12)) = 

q-®^ , (t u x u h)i; a2 (t2,x 2 ,i 2 )) , (C3) 

where B is the 2n x 2n matrix 



D B 
B l D 



(C4) 



with Z? given by (57) and 

B = 2crrl„ + a 2 $ t + t 2 1 



(C5) 



The eigenvalues the matrix D/4 provide the dimensions 
capturing the ipi — ip2 mixing. We will prove now that n 
of the eigenvalues of D/4 vanish and that the remaining 
n coincide precisely with the dimensions di given by (58). 
For this purpose we compute the characteristic polyno- 
mial det(P — x\2n)- First we move to the basis in which § 
has the form (45,46), performing the the transformation 











D B 
B l D 



D>* 

o* 



b<l 



D d B] 



(C6) 



In this basis D d is diagonal, whereas £?bd is block diago- 
nal. At this point we use the identity 54 



dct q^J = det(M) det(Q - PM^N) 



(C7) 



where M, N, P and Q are nxn blocks and M is invertiblc. 
Let us apply (C7) to det(D - xl 2n ) with D given by (C6). 
For Dd — xl n to be invertible we assume for the moment 
that x 7^ a 2 + t 2 + 2<TTSi with defined by (59). One 
gets 

det(D - a;I 2 „) = 

dct(L> d - xl n ) det{D d - xl n + B hd (D d - xI^BD . 

(C8) 

Being determinants of diagonal and of block diagonal ma- 
trices, the two factors in the right hand side of (C8) are 
easily computed. One finds 

n 

det(D d - xl n ) = ]J(x — (j 2 — t 2 — 2<jt Si ) , (C9) 



det(D d - xl n + B hd (D d - xln^BD = 
Y\?=ilx(x- 2n 2 -2t 2 -Aar Sl )] 

nr=i( a; - cr 2 - t 2 - 2oTSi) 



(CIO) 



Notice that the factor (C9) cancels precisely the denomi- 
nator of (CIO). Therefore, the characteristic polynomial 
we are looking for is 

n 

det(© - xhn) = Y[[x(x - 2o 2 - 2t 2 - AarSi)] , (Cll) 

which extends for any x by continuity and proves our 
statement. 



APPENDIX D: COMPOSITE TWO-FERMION 
OPERATORS 

We examine here the stability of the critical points 
under the perturbation with the composite operators 



$i(t,x,i) =: ip*ip 2 : {t,x,i) 



■ e iV^C-lv>i,R{vt~x)~ipi iL (vt+x)] . 



:, (Dl) 
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$ 2 (t,x,i) =: -02 "01 : {t,x,i) ~ 

. g-i\/^C-[ l Pi,fl(t't-2:)- l Pi,t(«t+a;)] . _ 

The relative two-point correlation functions are easily de- 
rived. One finds 

($t(h,x 1 ,i 1 )$i(t 2 ,x 2 ,i 2 )) = 
{$ 2 '{t 1 ,xi,ii)$ 2 (t 2 ,x 2 ,i 2 )) ~ 

[D(vt 12 - x 12 )f- Si ^ [D(vt 12 + x 12 )f- Si ^ 

[D(vt 12 - x 12 )]- C - S ^ [D(vt 12 + i u )]A^ , (D3) 

(<$>* 1 {t 1 ,x- L ,i- L )$ 2 (t 2 ,x 2 ,i 2 )) = 
($ 2 (t- L ,x 1 ,i 1 )$ 1 (t 2 ,x 2 ,i 2 )) ~ 

[X>(uti2 - £12)]^-^ [X>(uti2 + x 12 )]- c - 5i ^ 

p(vt 12 - £ 12 )] C " S ^ [X>(«ti 2 + xi 2 )f- S *^ ■ (D4) 

As before, the response of (D3, D4) under the scaling 
transformation (55) defines the 2n x 2n matrix 

D = C_ ("i 1 M ® (S + S* - 21) , (D5) 



the dimensions of the operators (Dl, D2) being the eigen- 
values of D/2. One easily finds that n of these eigenvalues 
vanish. The remaining n are given by 

* = C-(l-*i), (D6) 

where Sj arc defined by (59). Subtracting from (D6) the 
dimensions of the same operators on the line, one finds 
the nontrivial boundary dimensions 

^(boundary) = _ c 2 Sj _ (D?) 

At this point one can repeat the analysis performed in 
section VII for perturbations with a single fcrmion oper- 
ator. Comparing (62) and (D7) and using that > 0, 
we see that in the attractive regime g + < g_ the sta- 
bility properties of the critical points under the two dif- 
ferent perturbations are the same. In the repulsive case 
g + > g_ the behavior is inverted. The directions which 
were stable become unstable under perturbations with 
(Dl, D2) and vice versa. 
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